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A topological Dirac semimetal is a novel state of quantum matter which has recently attracted 
much attention as an apparent 3D version of graphene. In this paper, we report results on the 
electronic structure of the 3D Dirac semimetal Na 3 Bi at a surface that reveals its nontrivial ground- 
state. Our studies, for the first time, reveal that the two 3D Dirac cones go through a topological 
change in the constant energy contour as a function of the binding energy, featuring a Lifshitz point, 
which is missing in a strict 3D analog of graphene (in other words, Na3Bi is not a true analog of 
3D graphene). Our results identify the first example of a band saddle point singularity in 3D Dirac 
materials. This is in contrast to its 2D analogs such as graphene and the helical Dirac surface 
states of a topological insulator. The observation of multiple Dirac nodes in Na 3 Bi connecting via a 
Lifshitz point along its crystalline rotational axis away from the Kramers point serves as a decisive 
signature for the symmetry-protected nature of the Dirac semimetal’s topological bulk groundstate. 


Three-dimensional (3D) Dirac semimetals are materi¬ 
als whose bulk band structure features 3D Dirac fermion 
quasi-particles [H4l3| . It was theoretically predicted that 
a 3D Dirac semimetal state can be realized at the critical 
point of the topological phase transition between a band 
insulator and a 3D topological insulator in a space in¬ 
version symmetric bulk crystal [2] ■ The Dirac semimetal 
state realized in this way features a single 3D Dirac cone 
at the Kramers point where the bulk band inversion oc¬ 
curs [2). This type of 3D Dirac semimetal has been re¬ 
alized and observed in multiple material systems such 
as BiTl(S 0 . 5 Se 0 . 5)2 and (Bi 0 . 94 ln 0 . 06 ) 2 Se 3 l 13S, 15j] . How¬ 
ever, this approach requires fine-tuning the chemical dop¬ 
ing/alloying composition, which is difficult to control and 
also introduces chemical disorders limiting the mobility 
of the system. In 2012 and 2013, numerical band calcula¬ 
tions predicted 3D Dirac semimetal states in two stoichio¬ 
metric materials, Na 3 Bi and Cd 3 As 2 0 . (i . Unlike the 
fine-tuned case, the Dirac semimetal states in Na 3 Bi and 
Cd 3 As 2 show a pair of bulk Dirac nodes that are located 
along the rotation axis away from the Kramers point in 
calculation 0 0 - This new type of Dirac semimetal 
with multiple bulk Dirac nodes arises from the protec¬ 
tion by additional crystalline rotational symmetries 0 - 
[t] . They are believed to be stable within a range of 


material parameter space and robust to finite disorder 
or chemical doping. As a result, symmetry-protected 
Dirac semimetal states have soon attracted much in¬ 
terest Following these two band structure pre¬ 

dictions, several photoemission experiments studied the 
(001) surface electronic structure of Na 3 Bi and Cd 3 As 2 
16h 18I|. Although 3D Dirac cones have been observed 
in both materials 043 , it was found difficult to re¬ 
solve the existence of two bulk Dirac cones and their 
overlap/interplay in momentum space when studying 
the (001) surface. Therefore, the critically important 
unique features of the electronic structure of symmetry- 
protected Dirac semimetals have not been unambigu¬ 
ously probed in experiments. Revealing the true elec¬ 
tronic groundstate is fundamental for future works on 
any new quantum materials, either theoretically or ex¬ 
perimentally. This is even more important for a sym¬ 
metry protected Dirac semimetal such as Na 3 Bi because 
it is the phenomenon of multiple bulk Dirac nodes away 
from the Kramers point along the rotation axis that dis¬ 
tinguishes this new type of symmetry protected Dirac 
semimetal state from the former (fine-tuned) type that 
is achieved by fine-tuning (Tj. 

In this Letter, we solve this problem by studying the 
electronic structure of Na 3 Bi at a different surface - 
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(a) Na 3 Bi crystal structure (b) (100) crystal structure 


0 : Nal 0 : Na2 O Bi 


FIG. 1: Crystal structure and multiple Dirac nodes in 
NaaBi. (Characterization of NaaBi system, (a) Crys¬ 
tal structure of NasBi. The two Na sites and Bi atoms are 
marked with different colors, (b) Projected crystal structure 
at the (100) surface, (d) First-principles bulk band structure 
calculation, (d) Structure of bulk and surface Brillouin zone 
at (001) and (100) surfaces. Bulk Dirac nodes are marked by 
green crosses. Panels (c,d) are adapted from Ref. Q. 


the (100) surface by angle-resolved photoemission spec¬ 
troscopy (ARPES). The (001) and the (100) surfaces 
correspond to the (0001) and the (1010) surfaces in the 
hexagonal four-index scheme. Our studies, for the first 
time, reveal unambiguously that the electronic ground- 
state Dirac semimetal material NasBi consist of two 3D 
Dirac cones and that these two Dirac cones go through a 
Lifshitz transition, a topological change in the constant 
energy contour [19j| . which is missing in a strict 3D ana¬ 
log of graphene. Our results identify the first example 
of van Hove singularity in a 3D Dirac system, and there¬ 
fore pave the way for studying unusual phenomena re¬ 
lated to van Hove singularities in these novel materials. 
These observations are made possible because we were 
able to access the (100) surface. At the (001) surface 
(as the case in previous studies UlSl), the experiments 
were limited by the out-of-plane momentum resolution of 
ARPES probes. 

High quality single crystals of Na 3 Bi were grown by a 
method reported in Ref. [2l|. Ultraviolet spin-integrated 
ARPES measurements were performed at the beamline 14 
at the MAX-lab in Lund in Sweden and the ESPRESSO 
endstation installed at the beamline-9B of the Hiroshima 
Synchrotron Radiation Center (HiSOR) in Hiroshima in 
Japan. The energy and momentum resolution was bet¬ 
ter than 30 meV and 1% of the BZ for spin-integrated 
ARPES measurements at the beamline 14 at the MAX- 
lab and the ESPRESSO endstation at the beamline-9B 
of Hisor. Samples were cleaved in situ under a vacuum 
condition better than 1 x 10~ 10 Torr at all beamlines. 
Since bismuth trisodium is air sensitive, argon-filled glove 
boxes with residual oxygen and water level less than 1 
ppm were used in the entire preparation process. ARPES 
measurements were performed at liquid nitrogen temper¬ 
atures in the 14 beamline at the MAX-lab, and at 10-20 
K at the beamline 9B at HiSOR. Samples were found to 
be stable and without degradation for a typical measure- 
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FIG. 2: ARPES observation of multiple bulk Dirac 
cones, (a) ARPES constant energy maps obtained at the 
(100) surface using incident photon energy of 58 eV. Yellow 
dashed boxes indicate the surface Brillouin zone; notations of 
high symmetry points are shown at the Eb = 1000 meV panel; 
binding energies are marked in red on the lower left corner of 
each panel, (b) Selected ARPES k-E cuts, directions and 
rough fc-ranges of the cuts are indicated by black solid lines 
in the left panel. 


ment period of 24 hours. The first-principles calculations 
were based on the generalized gradient approximation 
(GGA) [12] using the full-potential projected augmented 
wave method as implemented in the VASP package. Ex¬ 
perimental lattice parameters were adopted from Ref. 
’21]. The electronic structure of bulk bismuth trisodium 
was calculated using a 9 x 5 x 5 Monkhorst-Pack fc-mesh 
over the BZ with the inclusion of spin-orbit coupling. 

NaaBi crystalizes in the hexagonal PQ^/mmc crys¬ 
tal structure with a = 5.448 A and c = 9.655 A j2Q] 
(a _1 /c _1 ~ 1.8). It is a layered structure consisting of 
a honeycomb atomic plane where the Na and Bi atoms 
are situated on the A and B sites respectively, sandwiched 
by layers of Na atoms which form a triangular lattice [see 
Figs. l(a,b)]. In the ionic limit, the bulk conduction and 
valence bands are expected to be the Na 3s and the Bi 
6p x ,y,z, respectively. However, first-principles calculation 
with spin-orbit interaction [5] shows that the conduction 
and valence bands are inverted at the bulk BZ center T 
point. Since the [001] axis (r — A direction in momen¬ 
tum space) is a C 3 rotational axis, all electronic bands 
along the T — A direction consist of states, which are 
eigenstates of the C 3 rotational operator. Interestingly, 
in NasBi, the lowest conduction and valence bands that 
go through the band inversion at the T point are found 
to have different C 3 eigenvalues @,0. This fact suggests 
that an energy gap cannot open between these two bands 
along the r — A direction even after including spin-orbit 
interaction. Consequently, two bulk Dirac band cross- 
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FIG. 3: Out-of-plane dispersion of the bulk Dirac fermions, a, ARPES constant energy maps along the k y -k z plane. 
Binding energies are marked on top of each panel. The k z [100] axis is built by measuring ARPES k-E maps with photon 
energies ranging from 40 to 76 eV. b, Schematic diagram of the k z mapping procedure. Since the cleaving surface is (100), k z 
corresponds to the [100] direction parallel to the crystallographic a axis, c, Linear dispersion along the third dimension (fc 2 ). 
Probing with different photon energies, the top of one of the bulk Dirac cones locate at different binding energies, constructing 
the linear k z dispersion, d, Linear k z dispersion seen from in-plane k-E cuts. It is clear from the data that the bulk Dirac 
cone does not touch the Fermi level except for hv = 58 eV, consistent with the three dimensional Dirac dispersion. For Panels 
c and d, top row presents raw data, Bottom row shows second derivative results along the momentum distribution curves. 


ings along the A — T — A direction on the opposite sides 
of the r point are shown in first-principles calculations 
[Fig. 1(c)]. Therefore, the two bulk Dirac points are ex¬ 
pected to be located along the [001] axi s [p , [fj, which 
is perpendicular to the (001) surface [16 h 18| . Since the 
out-of-plane momentum resolution of ARPES is much 
limited, it is difficult to resolve the existence of two bulk 
Dirac cones and their overlap/interplay at the (001) sur¬ 
face. By contrast, they project onto different points in 
the (100) surface BZ [Fig. 1(d)]. 

Fig. 2 shows our measured energy-momentum disper¬ 
sion of the two bulk Dirac cones at the (100) surface. 
The ARPES constant energy maps are shown in Fig. 
2(a). The surface BZ is found to be in a rectangular 
shape, which is consistent with the (100) surface termi¬ 
nation. The length ratio between the long and short edge 
of the surface BZ (yellow boxes) is measured to be about 
1.8 ~ a -1 /c -1 , also consistent with a (100) cleaving sur¬ 
face. Remarkably, at the Fermi level, our data clearly 
show that the Fermi surface consists of two unconnected 
band contours, which represent the two bulk Dirac cones. 
And the two bulk Dirac cones are observed to be equally 
spaced on the opposite sides of the T point along the 
A — F — A direction. A slight increase of binding energy 
up to 50 meV results in enlargement of the two Fermi 


pockets. At higher binding energies this two cones merge. 
The band contour further expands to form a neck region 
around the A points (top of the 3D BZ). The chemical po¬ 
tential seems to be slightly below the energy of the two 
Dirac points, since the observed two separated pockets 
still has some finite area. Fig. 2(b) presents several typ¬ 
ical ARPES k-E maps along various cutting directions 
[black solid lines in the left panel of Fig. 2(b)]. A clear 
linear dispersion is seen for the bulk Dirac cones in all of 
these maps, establishing the existence of massless bulk 
quasiparticles. In Fig. 3 we study the band dispersion 
along the out-of-plane direction that is perpendicular to 
the (100) cleavage surface, by varying the incident pho¬ 
ton energies of the synchrotron radiation. We expect the 
constant energy ( k y — k z ) maps to demonstrate an en¬ 
larging, circular-like band contour. From Fig. 3(a) one 
realizes that this is indeed the case. Moreover, the E-k z 
(S-fe[ioo]) cut shown in Fig. 3(c) is also consistent with 
the expected linear dispersion. In Fig. 3(d) the same fact 
is emphasized again by showing the k v -E maps for mul¬ 
tiple photon energies. One realizes from Fig. 3(d) that 
the Dirac band touches the Fermi level only at hv = 58 
eV. 

Since the two bulk Dirac cones are located very close to 
one another in momentum space, they inevitably touch 
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FIG. 4: Lifshitz transition and saddle point singularity 
in the bulk Dirac cones, (a) Lifshitz transition observed 
in ARPES constant energy maps. At Eb ~ 144 meV the two 
cones touch and hybridize to form a diamond shaped single 
band contour, marking the Lifshitz transition, (b) Theoreti¬ 
cal 3D constant energy maps at various binding energies, for 
comparison with Panel (a). Inset in the 134 meV panel shows 
band projections onto the (100) surface at the Lifshitz en¬ 
ergy. CEM: constant energy map. (c) Lifshitz transition of 
bulk Dirac fermions in NaaBi (top panel), compared with the 
Lifshitz transition of surface Dirac gas in a topological crys¬ 
talline insulator Pbo. 7 S 1 io. 3 Se (bottom panel), (d) ARPES 
k-E maps along Cut X and Cut Y defined in Panel (a) with 
and without the dashed lines. The dashed lines are guides 
to the eyes, which are drawn to only show the trend of the 
band dispersion at a qualitative level, (e) Saddle point sin¬ 
gularity observed from an anomalous increase of local density 
of states. The peak of integrated ARPES intensity marks the 
Lifshitz transition at Eb ~ 155 meV. The EDC is obtained by 
integrating the ARPES intensity at a (± 0.1 A 1 ,±0.1 A X ) 
window in the vicinity of the saddle point, which is the point 
where Cut X and Cut Y cross each other in panel (a). 


and hybridize with each other, giving rise to a topologi¬ 
cal change in the band contours, also known as a Lifshitz 
transition in the electronic structure. Fig. 4 demon¬ 
strates our systematic ARPES investigation on the Lif¬ 
shitz transition as well as the associated saddle point 
band structure. Fig. 4(a) presents a set of ARPES 


constant energy maps with higher momentum resolution 
than the one shown in Fig. 2. Besides the two bulk 
Dirac cones presented already in Fig. 2, we emphasize 
here the change of band structure where the two uncon¬ 
nected Dirac band contours merge. Experimentally, this 
Lifshitz transition happens at a binding energy of ~ 144 
meV, at which the two hole-like bands transform to a 
diamond-shaped loop. At higher binding energies, this 
loop enlarges while another concentric contour starts to 
appear at the zone center. This trend of band evolution is 
nicely reproduced in first principles band calculations, as 
shown in Fig. 4(b) where a set of 3D (/c[ 10 o], fc[oio], fc[ooi]) 
constant energy maps at selected binding energies are 
presented. The theoretical critical point of the Lifshitz 
transition is at£?B = 134 meV, only 10 meV away from 
the experimentally observed value. It can be seen in Fig. 
4(a) that the two cones touch at two special momenta, 
( k x , k y ) = (± 0 . 1 , 0 ) A .To directly identify the saddle 
point band structure as a result of the Lifshitz transition 
in experiments, we center our ARPES detector at one of 
the special momenta at ( k x , k v ) = (± 0 . 1 , 0 ) A and 
study two k-E maps, namely, Cuts X and Y as defined 
in the 3 rd panel in Fig. 4(a). As shown in Fig. 4(d), 
the hybridized band reaches its local energy maximum 
at ( k x , k y ) = (± 0 . 1 , 0 ) A along k x , while reaching a 
local energy minimum along k y . This behavior defines a 
saddle point singularity of band structure. At Fig. 4(e), 
an increase of local density of state (DOS) is observed 
by integrating the ARPES intensity around the saddle 
point, consistent with a local DOS maximum in a saddle 
point type of van Hove singularity in three-dimensions. 
In Fig. 4(c) we compare and contrast the Lifshitz tran¬ 
sition observed here with that observed in the surface 
states of a topological crystalline insulator Pbo. 7 Sno. 3 Se. 
In the n-typed sample of Pbo. 7 Sno. 3 Se, the two surface 
Dirac cones hybridize at Eb = 23 and 96 meV. The dif¬ 
ference between these two cases is that for NasBi, bulk 
Dirac quasiparticles, instead of the surface Dirac electron 
gas in the case of TCI, merge and hybridize. 


We note the following observations in our experiments: 
(i) we observe a pair of bulk Dirac nodes; (ii) they are 
located along the C 3 fc[ooil axis away from the T point; 
(iii) the distance between the two in momentum space 
is only 0.2 A ; (iv) the Lifshitz transition and saddle 
point singularity between the two are resolved. These 
observations are made possible because we studied ( 100 ) 
surface. At the (001) surface EMI , it is difficult due 
to the fact that the out-of-plane momentum resolution of 
ARPES is quite limited. We discuss the connection and 
distinction between the saddle point type of van Hove 
singularities in 2D and 3D band structures. Saddle point 
singularities in 2D Dirac electron gas have been exten¬ 
sively studied in twisted bi-layer graphene 23, 24] and 
the surfaces of topological crystalline insulators [25j. In 
2D, the DOS at the energy of the singularity diverges 
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logarithmically fashion leading to electronic instability 
and correlated physics. In contrast, in 3D, the DOS it¬ 
self reaches a finite maximum, but the derivative of the 
DOS does diverge. However, unconventional correlation 
physics can still arise if the bulk band structure is in a 
good quasi-2D approximation, such as the case in high- 
T c copper superconductors [26|. Despite this difference, 
our results have, for the first time, identify an example of 
band structure singularity in 3D Dirac materials. More¬ 
over, we note that strictly speaking, Lifshitz transitions 
are topological changes of the Fermi surface upon chang¬ 
ing some external parameter (doping, pressure, tempera¬ 
ture, etc). Thus to achieve such a Lifshitz transition and 
to access the singularity at the Fermi level, one needs 
to dope NaaBi to tune the chemical potential. A recent 
paper showed that controlling the chemical potential is 
indeed possible by finely adjusting the growth conditions 
of NaaBi [27]. Furthermore, we note that a very recent 
theory work has classified all possible 3D Dirac semimetal 
states Q: theoretically, there are three classes of Dirac 
semimetals Q, namely the (trivial) Dirac semimetal re¬ 
alized at the critical point of a topological phase tran¬ 
sition, the (topological) Dirac semimetal state with a 
pair of Dirac nodes protected by the rotational sym¬ 
metry, and the single-Dirac-cone Dirac semimetal state 
at the Kramers point protected by rotational symme¬ 
try. Our unambiguous observation of a pair of bulk Dirac 
nodes along the crystalline rotational axis away from the 
Kramers T point serve as a signature for the symmetry- 
protected nature of the Dirac semimetal state, and, for 
the first time, revealed that the groundstate of NasBi is 
a nontrivial 3D generalization of graphene. 
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